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Abstract 

Many-particle correlations due to Bose-Einstein interference are studied in ultrarelativistic heavy-ion 
collisions. We calculate the higher order correlation functions from the 2-particle correlation function by 
assuming that the source is emitting particles incoherently. In particular parametrizations of and relations 
between longitudinal, sidewards, outwards and invariant radii and corresponding momenta are discussed. 
The results are especially useful in low statistics measurements of higher order correlation functions. We 
evaluate the three-pion correlation function recently measured by NA44 and predict the 2-pion-2-kaon 
correlation function. Finally, many particle Coulomb corrections are discussed. 



I. INTRODUCTION 



Particle interferometry is an indispensable tool for the study of the space-time structure of the emitting 
source created in ultrarelativistic heavy-ion collisions. Recent experiments not only enable us to 

obtain a good spatial resolution in all 3 dimensions for the 2-pion correlation function, but also provide 
large enough data samples to investigate many-pion correlations Q . 

The basic issue one wants to address and resolve with three or more particle interferometry is, if there 
exist additional correlations beyond Bose-Einstein (BE) interference of identical particles from an incoherent 
source. In this paper we predict the three- and more particle correlation functions on these assumptions and 
explore whether the two-body correlation function suffices to describe the many-particle system. In particular 
we compare to recent data on the three-body correlation function. Any deviation found experimentally from 
this prediction signals new physics like, e.g., coherent emission or true higher order correlations. We use 
the measured two-particle correlation functions as input and therefore obtain a parameter free prediction of 
higher order correlation functions. 

Our starting point are the three-parameter fits to the two-pion correlation function Ci- This correlation 
function depends on the momentum difference, q, of the pion pair. In section 2 we provide a simple method 
for calculating many-body correlation functions due to Bose-Einstein interference between identical particles 
emitted incoherently from a source and discuss to which degree they can be determined from the two-body 
correlation function alone. To demonstrate our method, we reduce in section 3 the dependence on the three- 
relative momenta to a dependence on only the Lorentz invariant relative momentum and compare our result 
to the corresponding one-parameter experimental fit and data JjJ . The resulting expression is approximated 
by an expansion in the deviations of the pion source from a spherical shape. We show that the lowest order 
corrections are sufficient to reproduce the full expression. 

In section 4 we proceed to evaluate the three-pion correlation function C3 and compare it with recent 
results Q . We find that within the experimental errors the prediction based only on three-pion correlations 
can reproduce the data. Finally in section 5 we discuss Coulomb effects and present a simple way to correct 
many particle correlation functions. We conclude with a summary of our results. 

II. BE INTERFERENCE AND MANY-BODY CORRELATION FUNCTIONS 

We want to outline a general calculation of many-body correlation functions which can be done elegantly 
for relativistic heavy-ion collisions. The method reproduces in a simple calculation the two-, three- and 
four-body correlation functions found previously by diagrammatic techniques Jl2[ . Readers familiar with the 
results may skip this section. 

The approximations, that we will employ, are: 

• incoherent source, i.e., particles are emitted independent of each other and carry no phase information 
as is the case in (coherent) stimulated emission. This is a strong assumption, which can only be tested 
by comparison to experimental data. Existing two-body correlation functions are well described by 
fully incoherent pion and kaon emission (see, e.g., flio||). 

• plane wave propagation after emission. This assumes that Coulomb repulsion can be corrected for (see, 
e.g., [|l4| and section 5) and that other final state interactions are negligible. 

• particle momenta are much larger than relative momenta. This is a good approximation as the mo- 
menta of particles emerging from relativistic heavy ion collisions are larger than the typical transverse 
momenta which are ~ 300MeV/c. In comparison the interference occur for relative momenta, q, less 



• large multiplicity. The multiplicity in high energy nuclear collision is hundreds or thousands of particles 
which is much larger than the number of particles we want study correlations for. 



The N-body correlation function relates directly to measured cross-sections 

C(fcl ' fc2 '-'^ ) = ff d^ kl ....d^k N (1) 

where k\, k^ are the four- momenta of the N particles. 

The wave function for a particle emitted at the source point Xi with momentum kj is given by (xi\kj) = 
e lXik i . We used here the assumption that the particle propagates like a plane wave. We can then write down 
the symmetrized wavefunction for TV bosons, 

|V s (fci, .., k N )) = 7= E l^i* -> fc ^)> > ( 2 ) 

ViV. p 

where P is an operator representing all possible permutations of momenta insuring that the wave function 
is symmetric. Anti-symmctrization in case of fermions differs only by sign changes for odd permutations in 
equation (Jfy. The symmetrized wave function for 2 bosons is, for example, 

|V.(* lj fca)) = ^(|Ai,fe) + |*2 J fci)). (3) 

The symmetrized two-body wave-function in coordinate space is then 

(x 1 ,x 2 \ips(ki,k 2 )) = -L ( e lklXl e tk2X2 +e *^ e * 1X2 ) . (4) 
V2 

For an incoherent source, S(x,k), emitting particles at space-time point x with momentum k, the corre- 
lation function can be written as (see, e.g., Jll| for a derivation of the two-body case also with a partially 
coherent source) 

C N (ki,...k N ) = J S(xi,ki)...S(x N ,k N )\(xi, ...,x N \ip s (ki, ...,k N ))\ 2 d 4: x 1 d 4 x N , (5) 

when the source is normalized such that J S(x, k)d x = 1. 

Since ki ~ K = k\ + ... + k^)/N we can approximate S(xi, ki) ~ S(xi,K). We have furthermore 
permutation symmetry in the cc^ 's in Eq. (0). As a consequence we obtain for the wave function expectation 
value (for bosons) 

|(xi, ...,X N \l/) s (kl, fcAr))| 2 = — ^ (P'(h, — j k N )\xi, ...,X N ) (Xi, ...,Xjv|P(fcl, ...,k N )} 



p,p> 



_ g-txi(fci-P(fci)) e ix N {kN-P(k N )) ^ 

P 

We used here a symbolic notation. Every P(fcj) = fcj can be any of the particle momenta. However, the 
set P(fci), ...,P(fcjv) is the same as the set ki,...,kjf - just permuted around in all possible ways. For the 
two-particle case TV = 2 we obtain from Eq. (§) |(xi,x 2 |V's(fci, fe))! 2 = 1 + e^ 1 ^ 2 ^ 1 ^ 2 '. Note, that 
the permutation symmetry in coordinates Xi,...,xjsr insures that the permutations of either wave- function 
cancels the factor I/AH. 

From Eqs. (]|) and (||) we obtain 
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N 

C N (k u ...k N ) / S(x,K)e^- p(k ^d 4 x , (7) 

P i=l ^ 

where the sum is over all the possible permutations and P(fc,) = fej represents all possible momenta kj 
making up one permutation. Here, we can for few-body systems ignore factors of (n — N)/n due to the large 
multiplicity n in relativistic heavy ion collisions. 

It is now convenient to define the Fourier transforms 



F(q) = J dx S(x, K) exp(iqx) , (8) 

where the dependence on the total momentum K has been suppressed for later convenience. The TV-body 
correlation function is then 

JV 

c N (k u ...k N ) = j2Ii F ^> ( 9 ) 

P i=l 

where q^ = fcj — kj = ki — P(k,j) and the sum is again over all possible permutations. 

It is now straight forward to write down the first few correlation functions. The two-body 



C 2 (fci, ki) = 1 + F(q 12 )F(q 21 ) = 1 + |F(g 12 )| 2 ; (10) 



and the three-body 



C 3 (ki, k 2 ,k 3 ) = 1 + F(q 12 )F{q 21 ) + F(q 23 )F(q 32 ) + F(q 31 )F(q 13 ) 
+ F(q 12 )F(q 23 )F(q 31 ) + F (q 21 ) F (q 32 ) F (q 13 ) 
= 1 + \F(q 12 )\ 2 + \F(q 23 )\ 2 + \F(q 31 )\ 2 + 2Re[F(q 12 )F(q 23 )F(q 31 )} . (11) 

For the four-body correlation function we give only the final result 

C^fca.fcs,**) = 1 + \F(q 12 )\ 2 + \F(q 13 )\ 2 + \F(q u )\ 2 + \F(q 23 )\ 2 + \F{q 2i )\ 2 + \F(q 34 )\ 2 
+ 2 (Re[F( qi2 )F(q 23 )F(q 31 )} + Re[F(q 12 )F(q 24 )F(q 41 )] 

+ Re[F(q 13 )F(q 34 )F{q 41 )} + Re{F(q 23 )F(q 34 )F{q 42 )} ) 
+ 2 (Re[F(q 12 )F(q 23 )F(q 34 )F(q 41 )} + Re[F(q 12 )F{q 2A )F{q i3 )F{q 31 )] 

+ Re[F(q 13 )F(q 32 )F(q 24 )F{q 41 )} ) 
+ \F(q 12 )\ 2 \F(q 34 )\ 2 + \F(q 13 )\ 2 \F(q 24 )\ 2 + \F(q 14 )\ 2 \F(q 23 )\ 2 . (12) 

These expressions were also derived by a diagrammatic approach [ p"2| . 

In case of two-pion and two-kaon correlations the symmetrization/permutation should only be within each 
pair of identical particles and so the 2x2 correlation function becomes 

C 2 x2(q 7 r l7 r 2 ,q Kl K 2 ) = (1+ \F„ (q ffl W2 ) 1 2 ) (1 + |i^(q Kl K 2 )| 2 ) . (13) 

Since the pion and kaon sources may be different, F„ and Fx may also differ. 

The crucial assumption of an incoherent real source thus allows us to write down any N-body correlation 
function if we know the function F(q). Its norm can be extracted from the two-body correlation function but 
not its phase. If the Fourier transforms of the source distributions Eq. (||) are complex, F(qy) = |F(q;j)| e*^", 
the two-body terms |F(qij)| 2 in Eqs. (fol) and (O) are insensitive to the phase fcj, whereas the three-body 
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term -F(qi2)-P 1 (q23)-F 1 (q3i) in ( fi"l| ) generally will depend on the phase. However, translations in space and 
time by x leads to a phase change iqijX of -F'(qij) but due to momentum conservation, q\i + 523 + <Z3i = 0, the 
phase of F (a^y^F (a^)F (a^\) is unchanged. In general all correlation functions are invariant with respect to 
translations. Furthermore, if the source is symmetric, i.e., p(x) = p{—x), we find that -F(qij) are real from 
their definition Eq. (|J). If the source is not symmetric as, e.g.., when resonances contribute jfij, then the 
F(q;j)'s are generally complex and the phase in the last term of Eq. ( |TT| ) may carry interesting information. 
The phase can, however, not be studied by two-body correlations alone as seen from Eq. (|To|). 

III. TWO-PARTICLE BE CORRELATION FUNCTION 

In this section we analyse the two-body correlation function in relativistic nuclear collision and extract the 
crucial function for later use in many-body correlation functions. We will relate various parametriza- 

tions in particular as function of the invariant relative momentum as will be important for the three-body 
data. 

A commonly used parametrization of the two-particle correlation function is the three-parameter gaussian 
form 

C 2 (q h q s , qo) = l + A 2 cxp {-q 2 R 2 - q 2 R 2 - q 2 R 2 ), (14) 

where the components of the relative momentum q = pi — P2 are the longitudinal qi (along the beam axis), 
sidewards q s (perpendicular to both the beam axis and the pair momentum) and outwards q a (see Figure 1). 



El 
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q i3^k3 






J ■" ^ J z— longitudinal 


M y=sideward 



FIG. 1. View of the interaction region for the case of 3 particle emission. 

The dependence on K , or equivalently the rapidity y ~ y\ ~ j/2 or transverse momentum K± ~ p±\ ~ p±i, 
has been suppressed in Jl4|). It is implicit in the source size parameters, Ri tS ,o- Longitudinal expansion is 
revealed as a rapidity dependent Ri Jl~3|] while transverse expansion may be seen through R s<0 depending 
on the transverse momentum For an ideal chaotic source the parameter A2 would be unity but a 

number of effects like a partially coherent source, long lived resonances, final state interactions and Coulomb 
screening effects can lead to a smaller value for A2 pX( ] . Finally, we would like to point out, that it is common 
to include a term — 2q q\R 2 il for the parametrization in equation (|l4|). We did not include this term, since 
we concentrate in this paper on the experimental data of NA44. The generalisation of our calculations to 
include this term is straight forward. 

Due to low statistics it is commonplace to perform a gaussian fit of one variable only 
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C 2 (ginv) = 1 + A inv CXp ( - qf nv R f nv ) , 

where the invariant momentum squared is given by 



2 2 2 2,2,2 2 

9i„ v = q - 9t = <h + 9s + 9 - 4t • 



(15) 



(16) 



Here, qt = Ei — E 2 — ■ q/Ek since g <C If. Other variables, that also are used in parametrizations, are 
<Zfl =r = <Zs + 9o + <?? + 9? or = 9s + 9o + 9; 2 • I n the following we will perform an analysis in the variable 
qi nv but results can be taken over by changing the sign of (for q\ =T ) or removing (for q 2 ) the last term in 



equations ( |16[ ) and (|18|). 

Experimental data are often analysed in the longitudinal center of mass system (LCMS) of each pair, in 
which y = and K is along the outward direction. In this system q t = (3±q , where 



01 



yjm 2 + K 



(17) 



is the transverse velocity in the LCMS. The LCMS-frame and the different momenta are depicted in Figure 
1. The invariant relative momentum reduces to 



2 

9inv 



qf + q 2 s +q 2 o (1 -0±)- (18) 

We can relate C2(qinv) and 6*2(^1, q s ,Qo) if both are determined from the same data sample 



C 2 (qh 



J d 3 q C 2 (q\, q s , q Q ) S(qf nv - q{ 



9o 2 (l -01)) 



Jd 3 q6(qf nv -ql 



9o 2 (l-/3l)) 



(19) 



The (5-function enforces equation (|l6|), while the numerator assures proper normalization. 

In high energy nuclear collisions the sources sizes Ri, i=l,s,o, are often similar in size as in high energy pp, 
pA and S+Pb collisions (see [QJ^] or Tables 1 and 2). In such cases we can expand in the difference between 
the radii and obtain from (Tl9|) 



C 2 («h 



1 + A 2 exp (-qf nv R 2 



1 
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-iinv 



E 

—l.s.o' 



2\3\ 



{Ri~R z y + 0{qt v {Ri-R z ) 



where 



(20) 



R 2 = -(R 2 



Ri + R > 



is the average source size squared and 



R > =R /Jl~Pl 



(21) 



(22) 



is the outward radius boosted by the transverse velocity. In the limit of small (R 2 — Rf nv ) or small <7; nv 
we recover i?j nv = R from a comparison of (|20| ) and ([l5]). Notice that R is defined such that the leading 
correction in (|2tj| ) is second order in (R 2 — Rf nv ) and fourth order in <7j nv . 
Defining the transverse and longitudinal source deformation parameters 

p2 r?2 



2R 2 
R 2 



2R 2 



(23) 
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we can rewrite equation (|20|) as 

C 2 (<7inv) = 1 + A 2 cxp {-qf nv R 2 ) 



2 „4 d4 / 1 r:> 



1 + I 9 f nv i? 4 (- SI + S 2 ) + Oiq^R'S 3 ) 



(24) 



The definitions of these deformations or eccentricities are related to the deformation parameters (3 and 7 of 
Bohr & Mottelson §]. 

The deformations have a physical meaning. In many cylindrical symmetric models for sources with an 
emission time St, one finds that R 2 — R 2 = P 2 ± St 2 [||j^Jl^]- Thus 



St 2 = R 2 [S ± (2/3I 2 - 1) + Ji - 1]. 



(25) 



In the Bjorken 1-dimensional hydrodynamic expansion model, one finds to a good approximation that R 2 = 
Tt 2 I \J m 2 + K\ in LCMS, where r is the average invariant time it takes to freeze-out after the collisions 



and T is the temperature at freeze-out. In that case 

P± 



R 



(2Si + 1) 



(26) 



It may be a coincidence that the source appears to be sphericalB. Rephrased in terms of the deformations, 
we can say that the emission time is small and that the frcczc-out happens when longitudinal expansion is 
about as large as the transverse size of the system, i.e. when the source expands in three dimensions. 

We applied these results to the data samples taken in the experiment NA44 |Q,9. Figure 2 shows in the 
upper left corner the results for the reaction S + Pb — > tt + tt + X, while the upper right corner corresponds to 
the reaction p + Pb — > tt + tt + X . In the two plots in the bottom row of Figure 2 we performed the analogous 
calculations for kaons. The lower left shows the results for the reaction S + Pb — > K + K + X, while lower 
right corresponds to the reaction p + Pb — > K + K + X. A three-parameter fit was applied to the data samples 
with the results listed in Table 1 and Table 2. 





A 2 


Ri (fm) 


Rs (fm) 


Ro (fm) 


(px) (MeV) 


S+Pb 


0.56±0.02 


4.73±0.26 


4.15±0.27 


4.02±0.14 


150 


P+Pb 


0.41±0.02 


2.34±0.36 


2.00±0.25 


1.92±0.13 


150 



Table 1. Parameters of the three-parameter fits to the NA44 data sample at 200 AGcV for the two-pion 
correlation function (from H). 





A 2 


Ri (fm) 


Rs (fm) 


Ro (fm) 


(px) (MeV) 


S+Pb 


0.82±0.04 


3.02±0.20 


2.55±0.20 


2.77±0.12 


246 


P+Pb 


0.70±0.07 


2.40±0.30 


1.22±0.76 


1.53±0.17 


237 



Table 2. Parameters of the three-parameter fits to the NA44 data sample at 200 AGeV for the two-kaon 
correlation function (from B). 



*By spherical we mean Ri — R s = R' a but not necessarily a spherical and isotropic source 
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S+Pb— >7C 7U X 



P+Pb- >7C 7L X 




— equation (23) 

— experimental fit 
• data points 



— equation (23) 

— experimental fit 
• data points 




S+Pb->KK + X 



p+Pb->K + K + X 



equation (23) 
experimental fit 
' data points 





FIG. 2. Two-pion correlation function for the reaction S + Pb — > tt + -k + X (upper left), p + Pb — > 7r + 7r + X (upper 
right), 5 + Pfe -> (lower left) and p + Pb -> K + K + X (lower right). The data points are from NA44 [lJ|J, 

the solid line is calculated using (|24|) and the dashed line is the 1-parameter fit (Table 5 and 6). 

We would like to mention that to apply equation ([li]) and (^4j) to the three-parameter fits we have for 
some data samples to take into account the experimental cut-offs for low momenta. In some instances a 
bin size of Sq ~ 10 MeV/c was applied for each momentum projection and the corresponding two other 
momentum projections were summed over all momenta between and g cut ~ 20 MeV/c. This prescription 
modifies the value of A2. For the longitudinal direction we find, for example, 



r , \ f cut dq s /•*»* dq 
C 2 {qi) = / / C 2 {q\,q s ,qo) 

JO 9cut Jo Icut 



l + X 2 cxp(-qfRf) - 



2 2 7r Erf(i? s ? cut ) Erf(i? 9cut) 



4 Rs q C ut 



R q c 



1 + A 2 cxp(- (7l 2 i?i 2 ) (1 _ - i^i? 2 ), q cut Ri « 1, 



(27) 



where Erf(x) is the error function. Corresponding equations apply to the other two directions and result in 
an apparent reduction of the experimentally measured A2. In the case of NA44 this prescription was not 
applied in the fitting procedure ||. 

In Figure 2 we find reasonably good agreement between the data points (black dots), the 1-parameter fit 
(dashed line) and our calculation using equation ( p4| ) (solid line). The approximation (|24|) is in this case 
hardly distinguishable from the exact expression ([IJ]). The parameters used for the upper row of Figure 
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2 are listed in Table 3, while the parameters used for the lower row in Table 4. The parameters for the 
1-parameter fit are listed in Table 5 and 6 for pions and kaons respectively. 





#inv (fm) 


(fll) 


Si 


s± 


S+Pb 


4.98 


0.73 


-0.05 


0.35 


P+Pb 


2.41 


0.73 


-0.03 


0.34 



Table 3. Source parameters for the two-pion correlation function determined according to equations (|l 
and 





Rinv (fm) 


iPx) 


Si 


Sx 


S+Pb 


2.90 


0.44 


0.04 


0.18 


P+Pb 


1.84 


0.43 


0.35 


0.20 



Table 4. Source parameters for the two-kaon correlation function determined according to equations ( |l7j ) 
and (||). 







#inv (fm) 


S+Pb 


0.56 + 0.03 


5.00 + 0.22 


P+Pb 


0.38 + 0.03 


2.89 + 0.30 



Table 5. Parameters of the 1-parameter fit to the NA44 data sample for the two-pion correlation function 
(from |). 





Ainv 


-Rinv (fm) 


S+Pb 


0.92 + 0.08 


3.22 + 0.20 


P+Pb 


0.68 + 0.06 


1.71 + 0.17 



Table 6. Parameters of the 1-parameter fit to the NA44 data sample for the two-kaon correlation function 
(from |). 

IV. THREE-PARTICLE BE CORRELATION FUNCTIONS 

We can now apply the results of the previous sections to predict all three- and more particle correlation 
functions. The three-pion correlation function recently determined by NA44 Q is especially relevant. A 
1-parameter Gaussian fit was applied to the data and the radius and A-parameter of this fit are listed in 





A 


-Rinv (fm) 


S+Pb 


2.1 + 0.4 


2.7 + 0.2 



Table 7. Parameters of the 1-parameter fit to the NA44 data sample for the three-pion correlation function 
(from f§). 

Table 7. The data was fitted in terms of the sum, Q3, of the relative invariant momenta of the 3 particles 
labeled 1,2,3 
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Q2 2 i 2 i 2 

3 = Ql2 + 9l3 + 923 • 



(28) 



The relative invariant momentum of particle i and j is hereby given as in equation (|16[) (see also Figure 1) 



2 2 2 

<l <lr 9t,ij 



(29) 



(30) 



We proceed to evaluate the three-pion correlation function C 3 (Q 3 ), starting from 

QtQ\ = I dc * 12 ^23 C 3 (qi2, qi 3 , q23) 8(Qj ~ gig ~ gig ~ 923) <? 3 (qi2 + qi3 + Q23) 

/dqi2 rfqi3 c?q23 8(Q% - q 2 2 - q 2 3 - q\ 3 ) <5 3 (qi2 + q i3 + q23) 
The 5-function of the sum of the relative momentum vectors assures that the 3 particles span a triangle. 

As mentioned we cannot extract the phase from two-body correlation function. Wc shall therefore in the 
following assume that -F(qy) are real functions which then are known from the two-body correlation function. 
The three-body and any many-body correlation function is then completely given by equation ( |Tl| ) and (|^). 

At present statistics do not allow experimentalists to map out the three-body correlation function of all 
six variables (nine including K) and we therefore reduce it according to Eq. (j3(j|). 
In the case of a spherical source, we obtain from (]3C]) and ( p"l| ) 

(31) 



C 3 (Q 3 ) = 1 + 3 A 2 cxp{~x) h{x) + 2 A* 5 exp(-y) 



and for a slightly deformed source, we obtain 

C 3 {Q 3 ) = 1 + 3A 2 exp(-a:) 
+ 2 A* 5 exp(-y) 



/l(x) + l5 ^*+^ (4 /,(■<■>-*/:>(,■) 



1 + 7^(^1+^) 



0{R 6 Q 6 3 6^ 



(32) 



where we used the dimensionless variable x = R 2 Q§/3 and y = R 2 Q%/2. I v is related to the standard 
Bessel function l v 



i v { x ) = (~y i v ( X ) = £ 



(x/2) 2 " 
n\(n + v)\ 



Figure 3 depicts this result, together with a 1-parameter Gaussian fit of reference Qj. A satisfactory 
agreement is achieved. Based on the currently available data we seem not to need any correlations beyond 
the two-pion correlations. 
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Figure 3 



deformed s 

spherical s 

tit to data 



FIG. 3. Three-pion correlation function for the reaction S + Pb — > tt + tt + tt + X . The data points are from NA44 Q, 
the solid line is calculated using ( [3l| ) , assuming a spherical source (no deformations) and the dashed-dotted line uses 
the approximation (|3l|), allowing for deformations. Finally, the dashed line is the 1-parameter fit (Table 7). 

The three-body correlation function of equation (^Tj) or ( |32| ) is a sum of two gaussians. However, the 
current experimental resolution cannot resolve the last term in these expressions, since it decays with an 
3/2 times larger exponential than the 3 explicit two-particle contributions and the data does not explore 
sufficiently small Q3. More statistics is needed to resolve this second term which may be available in Pb+Pb 
collisions at the SPS or future RHIC or LHC heavy ion collisions. It is the second term that may carry 
a phase if, e.g., the source is asymmetric. It would therefore be most interesting to study three and more 
particle correlations at small Q3 to learn more about the source. 

The method can be generalized to study different aspects of the 9 dimensional phase space, if the ex- 
perimental resolution should permit so, or extended to calculate any many-particle correlation function of 
invariant momenta. 

As a last example, we would like to proceed to the 4-particle correlation function. Current data samples are 
not large enough to study pure 4-pion correlation functions but they do allow us to study for example two- 
pion-two-kaon correlation functions 6*2x2 of equation ([l3|). As described above the two-particle correlation 
studies for pions and kaons can be used to obtain the invariant radii and eccentricities for the two particles 
and subsequently employing equations (|24|) and ( ]l0| ) allows us to calculate 62x2(97™, 9kk)- If the data should 
show any deviation from this result, we would have evidence for cross-correlations between pions and kaons. 



V. COULOMB CORRECTIONS 



In a recent paper |14j | Baym and Braun-Munzinger show that the standard use of the Gamow correction 
for Coulomb repulsion between pions of like charge is incorrect. The relevant length scale for the two pion 
wavefunction is their classical turning point in their Coulomb field, ~ Aam^/q 2 ^ fm, for relative momenta 
<7^10 MeV/c, and not the two-pion Bohr radius, qq = = 387 fm. Since the turning point is smaller 
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than the typical source size in high energy nuclear collisions, the two-body interactions are insignificant as 
compared to those of the source. It is instead a better approximation simply to assume classical trajectories 
for the two charged particles in their Coulomb field, Vc(r) — ZiZ^e 1 '/r, once they are separated by a distance 
ro of order the source size. A good description of tttt and irp correlation data from nuclear collisions at AGS 
was obtained with an initial separation of ro ~ 9 fm. 

The simple model of Ref. ju| can be generalized to N particles. Energy conservation implies that the final 
energies, Ei, of the particles are related to the initial ones, by 

Y J E % = Y,E^ + '-— , (33) 

i—l i=X 

where Rq is the N-particle equivalent of r and corresponds to an average separation between the particles. 

Particle number and total momentum conservation relates the final N-particle distribution njv(ki, ■••jk/v) 
to the initial distribution n^(ko,i, ...,ko jv) by 

N N 



n N 



(k x , ...,k iV )rf 3 k 1 ... ( i 3 k i v5(iVK - k = ^(ko,i, .... k , JV )d 3 ko ll ...d 3 ko l3 5(2VK - £ k 0ii ) 



(34) 



In terms of relative momenta, qij = ki — kj, and average four-momentum, K — {k\ + ... + krj)/N, we find 
that the invariant momentum of N particles is 

JV N 

q^JE 4 = -( NK ) 2 ( N ™) 2 = (£ E *f ^ V 2 + k 2 ) > (35) 

i,j=l i=l 

and likewise for initial momenta. With energy conservation in the Coulomb potential, Eq. (|3^), we obtain 

Q 2 n-QIn = N 2 (N-1)E k ^-, (36) 



where Ek = TOj_coshY. Defining the corresponding N-body correlation functions analogous to ( |30| ) we 
obtain from Eq. @ and @ 

~. \ 3JV-5 



C N (Q N ) = I ^ ) C° n (Qo,n) ■ (37) 



The measured correlation function is thus the original one, C%(^jQ 2 N — N 2 (N — 1)Ek e 2 /Ro), scaled by 
powers of Qo } n/Qn- In the two-body case the results reduce to those in Ref. jL4j. As an numerical example, 
we consider the three-body correlations of NA44 in LCMS (where Y = 0) for pions with (mi) ~ 300 MeV. 
Assuming i?o ~ tq ~ 9 fm the three-body invariant momentum is Q 2 — Q 2 3 ~ (30 MeV) 2 . In comparison 
the lowest measured data point in [Q is Q 3 ~ 45 MeV. 

The data of NA44 [|]J|] was corrected by the Gamow factor 
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C2{qv2)=T{q l2 )C Q 2 {q 12 ) 

TV ^ 2nf l 
r (<7l2) = 77 s 7 

exp{2irri) — 1 



7] = e 2 m/qi2 . 



(38) 



(39) 



The data for the three-pion correlation function [Q was corrected ad hoc by a three-body Gamow factor 
r(gi2)r(<7i3)r(g23)- Hereby Coulomb effects are overcorrected for but so are they in the two-body case at 
the smallest measured relative momenta . A better analysis including Coulomb effects correctly in two- and 
three-body correlation functions should be performed with a realistic charge distribution and not just a one 



In this paper we develop a simple procedure to predict many-particle correlation functions from Bose- 
Einstein interference between identical particles emitted independently from an incoherent source. Using 
the best available fits for the two-particle correlation function, we calculate the higher order correlation 
functions. Any deviation from this result should signal interesting new physics. 

We apply our prescription to pion correlations. The experiment NA44 supplies a single data sample, which 
allows one to study two- and three-pion correlations, as well as two-pion-two-kaon correlations. We test 
our approach against the data for two- and three-pion correlations in terms of the total invariant momentum 
of the pions and find good agreement. In the case of three-pion correlations no additional three-particle 
correlations seem to be necessary to predict the available data. We show, how to evaluate the two-pion-two- 
kaon correlation function from the currently available data. 

Larger data samples, like the one from the current Pb+Pb run at CERN will provide us with larger data 
samples and allow us a more thorough investigation of this issue. This data sample will not only improve 
the resolution for the three-pion correlation function, but also allow us to study 4-pion correlations. 

Finally, Coulomb effects were discussed and a simple N-particle generalization of the two-particle Baym 
& Braun-Munzinger correction was given in terms of invariant momenta. 



We would like to thank Hans B0ggild, Rama Jayanti, Bengt Lorstad, Sanjeev U. Pandey and Janus 
Schmidt-S0rensen for hclpfull discussions and for providing and explaining the necessary data. 



parameter model as in Jl4| and Eq. (|3^). 



VI. CONCLUSION 
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